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The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, 
convected steadily past a horizontal flat plate and sedimenting under the action of 
gravity, is examined. It is shown that, in the limit as Re^oc and e >0, where Re is the 
bulk Reynolds number and e is the ratio of the particle radius a to the characteristic 
length scale L, the analysis for determining the particle concentration profile has 
several aspects in common with that of obtaining the temperature profile in forced- 
convection heat transfer from a wall to a fluid stream moving at high Reynolds and 
Prandtl numbers. Specifically, it is found that the particle concentration remains 
uniform throughout the 0{Re ~ li 2 ) thick Blasius boundary layer except for two 0(e 2/3 ) 
thin regions on either side of the plate, where the concentration profile becomes non- 
uniform owing to the presence of shear-induced particle diffusion which balances the 
particle flux due to convection and sedimentation. The system of equations within this 
concentration boundary layer admits a similarity solution near the leading edge of the 
plate, according to which the particle concentration along the top surface of the plate 
increases from its value in the free stream by an amount proportional to X 5/6 , with X 
measuring the distance along the plate, and decreases in a similar fashion along the 
underside. But, unlike the case of gravity settling on an inclined plate in the absence 
of a bulk flow at infinity considered earlier (Nir & Acrivos 1 990), here the concentration 
profile remains continuous everywhere. For values of X beyond the region near the 
leading edge, the particle concentration profile is obtained through the numerical 
solution of the relevant equations. It is found that, as predicted from the similarity 
solution, there exists a value of X at which the particle concentration along the top side 
of the plate attains its maximum value (j) m and that, beyond this point, a stagnant 
sediment layer will form that grows steadily in time. This critical value of X is 
computed as a function of 0 S , the particle volume fraction in the free stream. In 
contrast, but again in conformity with the similarity solution, for values of X 
sufficiently far removed from the leading edge along the underside of the plate, a 
particle-free region is predicted to form adjacent to the plate. This model, with minor 
modifications, can be used to describe particle migration in other shear flows, as, for 
example, in the case of crossflow microfiltration. 


1. Introduction 

We consider the steady laminar flow of a suspension of sedimenting monodisperse 
solid particles past a horizontal flat plate. Far from the plate, the free-stream velocity 
U w and particle concentration (j) s are both uniform. We also suppose that the particles 
are small enough for their Reynolds number, based on their diameter and either the 
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Stokes sedimentation velocity or the shear rate inside the viscous layer formed near the 
plate, to be very small. 

This flow situation constitutes a fairly general prototype for studying the effect of 
shear-induced particle migration in shear layers. It is also a natural extension of the 
problem considered by Nir & Acrivos (1990) in their study of sedimentation onto an 
inclined plate, in which a concentrated sediment flows along the inclined plate owing 
to gravity. In that case, the shear created within the flowing sediment gives rise to 
shear-induced particle diffusion which opposes the sedimentation flux and prevents the 
particles from accumulating onto the upper surface of the plate. This occurs within a 
thin region close to the wall, termed the viscous layer, where viscous forces are 
balanced by buoyancy. In fact, the structure of that solution has several points in 
common with that of free-convection heat transfer from a vertical plate at high 
Grashof and Prandtl numbers with the notable difference that the particle 
concentration profile is discontinuous across the interface separating the suspension 
from the flowing concentrated sediment. Nir & Acrivos (1990) derived expressions for 
the particle velocity and concentration profiles within this sediment layer on the basis 
of an effective continuum model, valid for very small particle sizes relative to the 
characteristic length scale of the flow, which incorporated the effects of shear-induced 
particle diffusion (Leighton & Acrivos 1986). Their model was recently extended by 
Kapoor & Acrivos (1995) who also tested experimentally a number of theoretical 
predictions pertaining to the thickness of the sediment layer and the corresponding 
particle velocity profiles, and confirmed the theoretical results which were arrived at 
through ab initio calculations that did not entail the use of adjustable parameters. 

In view of the reliability of the effective medium model for buoyancy-driven flow, it 
would be of interest to consider an analogous situation in which forced convection 
plays an important role. In the problem examined here particles sediment onto the 
plate, where they tend to accumulate forming a concentrated sediment. At the same 
time, owing to the interaction between the free stream and the plate, a Blasius-type thin 
shear layer is formed near the plate, which subsequently induces shear-induced particle 
resuspension and thereby reduces their degree of accumulation. In addition, this whole 
process is affected by the presence of convective particle transport. In fact, as will be 
seen, an estimate of the size of the region within which the particle volume fraction 
variation is confined can be obtained simply by balancing the corresponding convective 
and diffusive particle fluxes. 

The present study constitutes a first attempt towards constructing a more 
comprehensive model of particle migration in shear layers which could then be used to 
describe various processes of significant practical interest. Such an example is crossflow 
microfiltration, in which a suspension of neutrally buoyant rigid spherical particles is 
made to flow in a channel having a flat porous bottom surface, under conditions where 
a pressure drop is imposed between the suspension and the permeate side of the 
microporous membrane. Consequently, fluid is sucked through the membrane, 
whereas the particles accumulate along the surface, forming a thin concentrated 
particle layer. Moreover, if this suction velocity is small relative to that in the centre 
of the channel, then, for all practical purposes, the flow within the bulk of the 
suspension is fully developed and is not affected significantly by particle migration. As 
a result, a known tangential shear stress is exerted along the particle layer which, 
together with convection, will determine the concentration profile (Davis & Leighton 
1987). In a different context, the shear-induced diffusion of raindrops as they approach 
an airfoil can affect the structure of the thin water film on its surface whose existence 
may be responsible for the observed loss in efficiency of airfoils under heavy rain and 
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high wind shear (Dunham 1987). Detailed analysis of these two flow situations is left 
for a future investigation; however, the mathematical and physical similarities between 
them and the problem at hand are fairly evident. 

In §2 the problem is formulated mathematically on the basis of the effective-medium 
model, while in §3 the basic equations are reduced to their boundary-layer form. It will 
be seen that the structure of the mathematical system is similar to that which describes 
laminar heat transfer from a flat plate to a fluid stream moving past it at high Reynolds 
and Prandtl numbers. In §4, a similarity solution is constructed valid near the leading 
edge of the plate according to which the particle concentration within the diffusion 
layer approaches smoothly its constant value, in the adjacent momentum layer 
without experiencing the discontinuity found by Nir & Acrivos (1990) in the absence 
of forced convection. Finally, in §5, the results of a numerical solution to the 
boundary-layer equations are presented which extend the theory beyond the leading 
edge of the plate. It will be seen that the particle concentration along the top side of 
the plate increases monotonically from its value at the tip, f) s , until it reaches its 
maximum possible value, at which point a stagnant layer is formed the thickness 
of which will increase continuously in time. In contrast, along the underside of the 
plate, the particle concentration decreases monotonically until a particle-free region is 
attained. 

2. Problem formulation 

We wish to examine the flow of a suspension of non-colloidal spherical solid particles 
in a viscous fluid convected steadily towards a semi-infinite flat plate, as shown in figure 
1, and sedimenting under the action of gravity. The particles are supposed to be small 
enough that their inertia is negligible. In the following a will denote the particle radius, 
t/ x the free-stream velocity, and p s , and p s , respectively, the effective density and 
viscosity of the suspension in the free stream, where the particle volume fraction 
equals <j> s . The suspension will be modelled as an effective continuum Newtonian fluid 
with concentration-dependent physical properties. 

In this context the volumetric average velocity of the binary mixture is defined as 

u = u p <j) + u f {\ -0), (2.1) 

where u p , u f are the particle and fluid ensemble-averaged velocities at any location in 
the mixture. Owing to the incompressibility of the fluid and of the solid particles, the 
effective suspension density becomes 

p(<f>) = Pf+(j> P -pf)<f> = p s y(0X (2.2) 

where p f and p p are the pure fluid and solid particle densities respectively and y(^) is 
the effective density of the mixture rendered dimensionless with p s . Similarly, the 
effective viscosity p{<j>) is set equal to p s A(<fi). The equations governing the motion of 
the mixture are rendered dimensionless using L, as yet unspecified, as the characteristic 
length scale, [7 X = |Ly as the characteristic velocity, and the inertial scale p s U % for the 
pressure. Consequently, following Nir & Acrivos (1990) the continuity and momentum 
equations for the mixture become, respectively, 

= 0, (2.3) 

and Rey(<fi)u-Vu = — Re VP + V •[A((f)(Vu + Vu T )] + -t l y^(<f> — 0 S ), (2.4) 

gIRe L 

where g is the gravitational acceleration. Besides the dimensionless parameters 
characterizing the flow are: the Reynolds number. Re = U. x p s L/p s , the Grashof 
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Figure 1 . Schematic of the flow. 


number, Gr = \{p p — p f )gp s L? / fi], and the size ratio e = a/L\ e 1 for the effective- 
medium model to be valid. The hydrostatic pressure, (p p — p f ) <fi s gy, has been 
incorporated into the pressure field P , which now vanishes far from the plate 
(henceforth, x and y will indicate, respectively, the Cartesian coordinates along and 
perpendicular to the plate as shown in figure 1). 

In dimensional form, the total particle flux is given by 


/ = -£>, -- |V«|) + (i i* + u) <j>, (2.5) 

with D J = |Vw| a 2 /3(<j>), D 2 = \ Vu\ a 2 K c <j> 2 , (2.6) 

where |Vw| = {Vw: Vw} 1/2 is the norm of the local gradient of the bulk velocity field 
(Leighton & Acrivos, 1987a, £), «* is the particle slip velocity, 




9 n f A(<p) 


(2-7) 


and ji f is the pure fluid viscosity. The first two terms on the right-hand side in (2.5) 
account for the diffusive flux of the particles due to concentration and shear rate 
gradients, respectively (Leighton & Acrivos 1987 a, b), whereas the last term 
corresponds to the convective flux. In addition, the particle slip velocity, given by the 
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first equality in (2.7), is set equal to the Stokes terminal velocity u t multiplied by a 
hindrance function /($$) which accounts for the influence of particle interactions. The 
expression for fl<j>) shown in (2.7) is that used by Nir & Acrivos (1990). The expressions 
given above for D ly D 2 and w* presuppose that the particle Reynolds number, whether 
in terms of the sedimentation velocity or the local bulk shear rate, is vanishingly small. 

In terms of the dimensionless variables introduced earlier, the steady-state particle 
balance equation becomes 


U -V(j) = e 2 V 


m i v «i n + vow r m- % m) 


( 2 . 8 ) 


At the solid boundary, the no-slip condition is applied together with the requirement 
that the total particle flux at the plate be equal to zero, i.e. 


u = 0 and 


y?(0)|Vw| V$4 + 




V(A(0)|V«|) 


gGr 

g Re 


m) 


■n = 0, 


(2.9) 


where n is the unit vector normal to the surface. At large distances from the wall the 
flow is uniform and the suspension is well mixed, hence 

u=U„, (2.10) 

In performing the calculations to be reported below, the relative effective viscosity 
of the mixture was taken to be of the form 


m 


/'W 

/<* 


(1 + 1 - 50/1 zllM. 
(1 + 1 . 5^/1 -MK)* 


4> m = 0.58 


( 2 . 11 ) 


(Leighton & Acrivos 1987 ft), where <j> m is the maximum possible particle volume 
fraction in the flowing suspension. Also, following Leighton & Acrivos (1987 ft), we let 

ftp) = i0 2 (l+ie 3 * * * * 8 ' 8?i ) and K c = 0.6. (2.12) 


3. Boundary-layer structure when Re §> 1 and e <? 1 

We consider the case of uniform flow past a flat plate at zero angle of attack. For 

large values of Re, based on the free-stream velocity U x and the length scale L, an 

0(Re U2 ) vorticity layer is formed near the plate. Above this region, the flow is uniform 
and the suspension is well mixed. Applying the standard Blasius-type boundary-layer 
scaling transformations 

U = u , V=vRe l/2 , X=x , V = yRe m , (3.1) 

with u and v denoting the longitudinal and transverse velocity components, respectively, 
we obtain, for the continuity, momentum and mass transfer equations, 




£(/ dV n 

1 — 0 , 

ex dY 


( W cU\ fl 


( U M +v w)-e? 



dP 

'ex' 


dP 
0 Y 


9 Grtf-P) 
2 Re h ' 2 ’ 


(3.2) 

(3.3) 

(3.4) 


U^+V^- 

ex of 




W I 


0 Y 


A(<f>) 0 Y 


H<t>) 


eu\ 


0F 


+ 


e 2 Gr i 
Re m ?-Y 


(3.5) 
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as well as the zero-flux condition at the plate, 


y= o, m 


W I 


8F 


# , Kctf 5 

^ T I" 


ay a(^) ar 


m 


at/| 


ay 


+*M)j? = o. 


Next, we examine the magnitude of e 2 Re 312 which can be rearranged into 


e 2 Re 312 = ( y 

[L 


a. 


W f2 = 


ui /2 


Ps / 


L 1/2 i>] /2 1 


UJLUJv s ) lf2 la 2 

L \v s 


(3.6) 


(3.7) 


where v s = p s /p s . Now, the term in square brackets is a measure of the shear rate near 
the plate as a result of the existence of a Blasius-type boundary layer, hence the entire 
expression corresponds to the shear particle Reynolds number within the boundary 
layer. As mentioned in the previous section, this was already assumed to be a small 
number in the context of our analysis. Having thus established that e 2 Re 212 is very 
small whenever the particle Reynolds number is sufficiently small, we can proceed with 
the analysis in a fashion similar to the case of heat transfer to a fluid stream convected 
past a flat plate at high Re and Pr . 

As e 2 Re w -^ 0, particle diffusion, both below and above the plate, is confined within 
a narrow layer attached to the plate which is embedded within the momentum layer as 
0. But, from the Blasius solution of the momentum boundary-layer equations, we 
know that, as 7^0, U ~ F, V ~ Y 2 (Schlichting 1968), and therefore, by requiring that 
the diffusive terms in (3.5) balance those due to convection, we conclude that the 
thickness of the diffusion layer must be 0(e 2/3 Re 1/2 ). This in turn leads to the 
transformations 

Y ~ U V 

Y= ^ V = (3.8) 


e 2f3 Re lt2 ' 


e 2/3 Re l/2 ' 


(e 2/ W /2 ) 2 ’ 


so that within the diffusion layer we obtain, respectively, for the momentum and 
continuity equations 



a 

ay 

[ A( ^J 

= 0{e 2 Re 3 ' 2 ), 

(3.9) 



K-o\ 

ay 

(-£■)• 

(3.10) 

and 


W dV n 
xpd — ~ — 0. 
dX 0y 

(3.11) 


In the next section we will see that, in the context of our analysis, e 2l3 Gr/Re 2 ^ 0, 
which results in a uniform pressure distribution. Further, in view of (3.9), the shear 
stress remains constant across the diffusion layer and hence shear-induced particle 
diffusion takes place only because of gradients in the particle concentration. In terms 
then of the new coordinates, the particle balance equation becomes 


0 

3* ay eyL 


st/| 

ay 


d</> 


e 2l3 Gr 0 
h — .. -wo 


ayj Re m 0 y 


(3.12) 


with boundary conditions 


y = o, 


m 

Y- 


0 U 


ay 

+ 00, 


T Y Re 3 ' 2 




(3.13) 

(3.14) 
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where the negative F-values denote the region underneath the plate. We note that, in 
contrast to the conventional cases involving mass transfer across boundary layers, here 
the diffusion coefficient, which appears both in the particle balance equation (3.12) as 
well as in the boundary condition (3.14), is proportional to the local shear rate rather 
than being constant. This, as we shall see presently, affects the solution of (3.12) in a 
significant way. 


4. Similarity solution valid near the leading edge of the plate 

As e 2l3 Gr/ Re 312 0 particle transfer due to sedimentation vanishes so that 

and the Blasius velocity profile applies throughout the boundary layer. Therefore, we 
construct the leading-edge expansion 


0 — 0* + 


e 2l3 Gr 

Re 312 


W{X, P) + ... 


(4.1) 


which, when substituted into (3.9)— (3. 1 2), yields, to 0(e 2l3 Gr/ Re 312 ), 



(4.2) 


W csV n 
_+_ = 0, 
dx er 

-cs*? ~6¥ / e 

U^+ V-Z = 
dx ey H Ys 0* 


' sc/ 

dF 

. a? 

0 Y_ 


0£7 

Y = 0, /?(&)--, 

?-*±oo, IF -*•<). 


0f 


+<f>M = 0 , 


(4.3) 

(4.4) 


(4.5) 


We see then that, to this order of approximation, the shear rate within the diffusion 
layer, dU/dY, is everywhere equal to ±0.332 X~ 1/2 , as given, for y->0, by the Blasius 
solution for the corresponding constant-property case. In addition, the term containing 
the settling rate only appears in the boundary condition at the wall. It is this constant 
term in (4.5) that accounts for the increase in the particle concentration within the 
diffusion layer, or for the particle depletion below the plate. 

Letting, then, 


U 


, 0-332 9 

— Vl/2 


and V = + 


0.332 Y 2 
4X W ’ 


Y> 0 


(4.6) 


in (4.2)-(4.5), we reduce the problem to a PDE whose form along with its boundary 
conditions suggests the existence of a similarity solution. Indeed, by introducing the 
transformations 


0.332 [fl(<f>s)] m 


(4.7) 


we obtain the ODE 




12 dr} 


d 2 J\ 

d )f ’ 


V ^ 


(4.8) 
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Figure 2. Similarity solution obtained by solving (4.8) for positive rj. 


with boundary conditions 

¥ = -l as (4.9) 

as 7j->± oo, (4.10) 

which can readily be solved numerically to give F^rj), shown in figure 2. In addition it 
can easily be seen that F x {rj) = — F x { — tj) with a discontinuity at rj = 0. Also /^(0 + ) = 
1.41. 

To a first approximation then, within the diffusion layer, whose thickness relative to 
the length scale L is 0{Re~ xl2 e w Re x/2 ) = 0(e 2/3 ), the particle concentration is given by 

< 4II > 

implying that, along the top face of the plate, the increase in <f> from its value in the free 
stream is proportional to X 5/ *, whereas, below the plate, <f> decreases in a similar 
fashion. It should also be noted that, unlike the similarity solution given by Nir & 
Acrivos (1990) for sedimentation over an inclined plate, here the particle concentration 
profile inside the diffusion layer, both above and below the plate, approaches 
continuously (f > s , its value in the undisturbed flow. In fact, an asymptotic analysis of 
(4.8) reveals that, as oo, F x {ij) ~ e^’^ 36 . The reason for this difference in the two 
concentration profiles, is that, in our case, the shear rate remains constant at any given 
X throughout the diffusion layer, whereas in Nir & Acrivos (1990) the corresponding 
shear rate is a monotonically decreasing function of Y and vanishes at infinity. 

Up to this point, the characteristic length scale L used in our analysis has been left 
unspecified. In what follows it proves convenient, however, to choose L such that 

c 2/z Gr &/(&) 

Re w 0.332 [A&)] 2/3 


(4.12) 
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giving 


~ 0.332[/?(&r 3 (aUV* rip., 

_ &/(&) V v» ) a 3 gAp 


in which case (4.11) becomes 

</> = 4> 8 + F 1 (v)X 5/6 + ... as x->o. 
We note that if L is chosen according to (4.13), 


a 3 g Ap\ 3/5 / J7 X a^ 3/5 


(4.13) 

(4.14) 

(4.15) 


which, for small particles, is normally very small, as required by our analysis. In the 
context of crossflow microfiltration Davis & Sherwood (1990) obtained a similarity 
solution that has some of the characteristics of the one presented here. For instance, 
in their analysis, the thickness of the sheared concentration boundary layer is also 
proportional to X 113 . Of course, their study is different in that they neglect gravity 
settling, they consider fully developed flow in the bulk of the suspension as opposed to 
the Blasius-type vorticity layer considered here, and finally they focus on that portion 
of the membrane where the particles form a thin stagnant sediment layer, cake, with 
concentration (j> m = 0.58. However, if one focuses on the entrance portion of the 
membrane where the flow is not fully developed and the cake layer has not yet formed, 
<f) w < <p m , and provided the transmembrane pressure drop is much larger than the 
down-channel pressure drop, in which case the wall permeation velocity, v w9 is taken 
to be constant, it can be seen that, as A^O, a similarity solution exists which is exactly 
the same as the one presented here with L defined such that 


Z h 

e* ,3 Re in O.332[/?(0 S )] 2/3 


(4.16) 


and a^^^Re- 112 , e 2 Re 312 < 1. (4.17) 

Substituting in (4.16) the parameter values used by Davis & Birdsell (1987), 
<f) s = 0.02, t/ x , = 3.7 cm s" 1 , v w — 3.0 x 10 -3 cm s -1 , a = 2 x 10 -2 cm, // = 20 cP, p — 
1.18 g cm -3 , in their experimental investigation, we obtain L « 2 cm, which conforms 
to the conditions in (4.17). When X ~ L the similarity solution is not valid and the 
full equations have to be solved; they are only slightly different from the full equations 
derived in the next section to provide the location X cr for which the stagnant cake layer 
forms. This is a very important design parameter as it controls the plugging of the 
channel. 


5. Solution valid when X ~ 0(1) 

Introducing the further transformation 

Y = W (&)] 1/3 

we obtain for the momentum equation (3.9) 



(5.1) 

(5.2) 
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which, when combined with the asymptotic form of the Blasius solution as 
by (4.6), leads to 


u=±m s )} 


1/3 


0.332 r d t 


* 0 


* i/2 Jo m' 


Y^O given 


(5.3) 


and 




(5.4) 


on account of the continuity equation. Substituting the above into (3.1 1)-(3.13), we 
then obtain for the particle concentration equation 


1 [ C Y 
± ^Uo 


d? 




w+. 1 


dX~2X 312 


r d 

Jo m l 


2A'(dA/d < 6)(0^/0X)]) 0^ 

+ Jj0F 


0 

0f 


m 


X ll2 dY 


) , 


+ Fm, y^O, (5.5) 


valid for all X , subject to the boundary conditions 


with 


<j> s as f oo, as X^O, 


4>M.) 


A&)A(0- 


(5.6) 


(5.7) 


As A^O the similarity solution can easily be recovered. It is also not difficult to show 
that the first correction to the similarity solution will be of the form F 2 (tj) A" 10/6 . 

Equation (5.5) has to be solved numerically. Moreover, in view of the form of the 
similarity solution as X^ 0, it is useful to rewrite equation (5.5) in terms of rj = Y/X 1/z 
and X , and, in order to avoid large longitudinal derivatives of the concentration near 
X — 0, to introduce 

H X (5.8) 

as the dependent variable. The resulting equation is parabolic in X and was discretized 
along then ^-direction using the Galerkin finite element method with the B-cubic 
splines as basis functions. Integration by parts eliminates second-order derivatives and 
increases the accuracy of the method to 0(h z ) where h is the size of the largest element. 
The longitudinal direction A" is a time-like coordinate and the integration is performed 
implicitly via the second-order trapezoidal rule. This gives rise to a set of nonlinear 
algebraic equations that are solved using Newton’s iterations. Unfortunately, because 
00/0A" appears inside the integrand in the second term of (5.5), the system matrix is 
fully populated and the CPU time increases like A 2 , where N is the number of 
unknowns. Typically 4-5 Newton’s iterations were sufficient to reduce the error, 
defined in terms of the Euclidean norm | Y (k+1) — Y k | 2 , to less than 10“ 8 ; Y k denotes the 
fcth iterant of the vector of unknown coefficient of the B-spline representation of the 
dependent variable H. The computation was terminated when the particle volume 
fraction reached its maximum value, <p m , on the top face of the plate, or when it 
vanished at the bottom. The results presented below were checked with respect to mesh 
refinement and are accurate up to, at least, the third significant digit. 



Forced convection and sedimentation past a flat plate 


311 



Figure 3. Comparison of the evolution of <j > w , the particle concentration along the top face of the 
plate, for (f> s = 0.2, as obtained via similarity {<j> s + 1.41T 5 ' 6 ) and as calculated via the exact solution. 
The computation stops when (j> w reaches its maximum value 

5.1. Results and discussion for the region above the plate 

The numerical solution of (5.5) and (5.6), to be referred to henceforth as the exact 
solution, generates a concentration profile, the shape of which is qualitatively similar 
to that given by the similarity solution. In addition, <j ) w , the particle concentration 
along the top side of the plate, is found to increase monotonically from the leading edge 
and to become equal to at a value of X which, according to our model, is only a 
function of <f> & . Beyond this point a steady-state solution of the system of equations 
given above does not exist and a stagnant sediment layer, having a uniform particle 
volume fraction <j) m , is predicted to form whose width continuously increases in time. 

The evolution of <fi w , the concentration at the top face of the plate, both as predicted 
by the similarity expression, 0,+ 1.4Uf 6/e , and from the exact solution, is shown in 
figure 3 for <j> s = 0.2. It is seen that the two results agree only up to X « 0.05, thereby 
implying that the range in X within which the similarity solution applies is very limited. 
This interval can be extended by adding to the result of the similarity solution the first- 
order correction, F 2 (rj = 0) X 101 *, but even so, with increasing distance from the leading 
edge, the discrepancy between the two sets becomes rapidly more pronounced. Indeed, 
as shown in figure 3, the similarity solution significantly overpredicts the rate of particle 
accumulation on the plate. In addition, in view of (4.1), (4.7) and (4.12), we can clearly 
see that the range in X over which the leading-edge solution applies decreases further 
when <f> 8 is reduced. 

On comparing the similarity solution with the exact solution two important factors 
are worth noting. Firstly, the shear-induced diffusion coefficient, /$(</>% is larger than 
which is the value used in the similarity formulation, and this reduces the rate 
of particle accumulation on the plate compared to that predicted by the similarity 
solution. On the other hand, the relative viscosity within the suspension, A(0), is larger 
than A(f s ) in the exact formulation and consequently the shear rate is smaller, which 
would lead to a higher rate of particle accumulation along the plate. As will be seen 
later, however, the effect of the increase in the shear-induced diffusion coefficient 
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Y Y 

Figure 4. Evolution of the particle concentration profile in the region above the plate, as calculated 
via the exact solution at different Allocations downstream, (a) <p s = 0.2 and X = 0.4, 0.5, 0.64, 0.78, 
0.92; X cr = 0.93 for this value of <j> 8 . ( b ) 0 g - 0.05 and X = 0.5, 1.1, 1.7, 2.5, 3.0, 3.5, 4.0, 4.41, 4.49; 
X cr = 4.5 for this value of 


dominates particle migration near the plate, hence the exact solution leads to smaller 
values for cj > w . 

The above argument loses its validity as the wall concentration approaches the 
maximum value 0 m , which is a singular limit in the sense that a very sharp increase in 
the concentration gradient, d<ft u ,/dX^oo, is observed. This is a result of the rapid 
increase of the viscosity; A(<^)^oo as Consequently, a region of progressively 

smaller shear rates is formed and particles continue to accumulate on the plate as they 
keep sedimenting (figure 3). Once the concentration reaches its maximum value, a 
stagnant sediment layer is predicted to form which will remain attached to the plate 
and grow in time monotonically. In their study of sedimentation over an inclined plate, 
Nir & Acrivos (1990) reached a similar conclusion when the angle of inclination fell 
below a certain threshold value. In the present case then, there exists a critical value of 
the plate length, X cr , decreasing in magnitude with increasing <j> s , beyond which a 
steady-state solution of the problem does not exist. It is understood that, in the vicinity 
of the region where <p w = </> m , the presence of the stagnant layer alters the geometry of 
the problem and hence the boundary-layer analysis presented in this study will no 
longer apply. 

Figure 4(a) shows, for <f> s = 0.2, the evolution of the concentration profile in the 
region above the plate as predicted by the exact solution. The most important new 
feature of this exact solutions concerns the appearance of an inflection point in the 
concentration profile, which is seen to be moving away from the plate with increasing 
distance from the leading edge. The reason for this is as follows. By evaluating (5.5) at 
the plate we obtain 

^ ^)+fw|] = ° at F = °> (5-9) 

which, when combined with the boundary condition (5.6), gives 


a 

af 


'80 

0 Y 



= 0 


0 2 0 _ F A(F/A) 
dY 2 A d0 


at 


r = o. 


(5.10) 
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But 


m 

m 


(l-^)&(2 + e 8 - 8 *«) 
( l -^)^(2 + e 8 ^) 


with 


$ = 4>w, 


(5.11) 


which decreases with increasing wall concentration, owing to the rapid increase in the 
shear-induced diffusion coefficient with \jr. Consequently, 0 2 ^/0F 2 < 0 at the plate. But 
since the concentration gradients are everywhere negative and become zero at the edge 
of the diffusion layer, the second derivative of the concentration profile must become 
positive at some distance from the plate, thereby explaining the existence of the 
inflection point. 

Referring to figure 4(a) again, it should be noted that, with increasing X , the 
inflection point moves further away from the plate and that the thickness of the 
diffusion layer, i.e. the distance along the F-direction over which the particle volume 
fraction decreases from <j> w to <f> 9 , increases faster with X than predicted by the similarity 
solution. In addition, the concentration profile can be viewed as consisting of two 
parts. Within the first, extending between the plate and the inflection point, the 
concentration gradient is relatively small near the plate but increases significantly as 
the inflection point is approached. It will be seen later that in this region the 
sedimentation flux balances shear-induced resuspension, while convection remains 
negligibly small. Consequently, the variation of the particle volume fraction is 
determined primarily by the magnitude of the shear-induced diffusion coefficient. In 
the second region, extending from the inflection point up to the edge of the diffusion 
layer, the concentration gradients decay to zero exponentially fast as far from 

the plate. 

At this point it will be of interest to discuss the predictions of our model as the free- 
stream particle volume fraction, <f> s , varies between values close to zero and <p m . The 
evolution of the concentration profile, as computed from the exact solution for cj) s = 
0.05, is shown in figure 4(b ) . On comparing this figure with figure 4(a) which applies 
when <f> s — 0.2, one can see that the concentration profile remains qualitatively 
unaffected by the variation in On the other hand X cr , the length needed for <j) w to 
equal <j> m on the top face of the plate, is significantly larger and the scale in the F- 
direction, over which changes in the concentration take place, is also larger. At the 
same time the concentration gradients far from the plate are larger than those observed 
in figure 4(a). In view of these two effects, and keeping in mind that the numerical 
discretization used in both cases was the same, the numerical error as X-+ X cr and in 
the vicinity of the inflection point, where the gradients are the steepest, is very large 
when <f> s = 0.05. In fact, it was found that, as 0 9 ->O the quality of the numerical 
solution in this region deteriorated, for the reasons cited above, so that in order to 
improve the accuracy around the inflection point it would have been necessary to 
introduce some form of regriding. This was not pursued in this study, since the 
accuracy of the numerical solution near the wall, where the concentration gradients are 
small, was not affected. Hence, we were able to obtain reliable results for the variation 
of the particle volume fraction at the wall, 0 MJ , for <f> s as small as 0.01. 

In view of the difficulties encountered in trying to compute a numerical solution as 
(j) s -> 0, it seems advisable to investigate the asymptotic structure of the solution of the 
exact equations in this limit. Close examination of figures 4(a) and 4(6) reveals that the 
variation of the particle volume fraction in the transverse direction, F, is determined 
primarily through the balance between sedimentation and diffusion in the region near 
the wall. As was already mentioned in the previous paragraph but one, the 
concentration gradient, 00/0 F, is small near the wall where the shear-induced diffusion 
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coefficient is large, whereas far from the wall 8^/6 f increases significantly as the 
particle volume fraction, and consequently the shear-induced diffusion coefficient, 
decreases (see also (5.11)). This so-called inner region occupies almost the entire 
sediment layer and determines its thickness. 

In this context, and guided by the behaviour of the numerical solution for small <fi s , 
we can recover, to leading order, the asymptotic structure of the exact solution as 
(j) s 0. More specifically, in this limit, we rescale the coordinates within the inner 
region according to 

X = X<j>*, Y=Yfo a 9 0> 0, (5.12) 

and express the particle volume fraction, which is an 0(1) quantity, by means of 

<f> = &o(X, Y) + <j> 8 <P x (X, ?HO(& 2 ). (5.13) 

We next make use of the exact integral particle balance, as obtained by integrating 
(5.5) across the diffusion layer and from X = 0 to X, 

< 514 > 

plus the fact, guided by the numerical solutions, that within the inner region the 
diffusion and sedimentation terms in (5.5) are of comparable magnitude and dominate 
those due to convection as 0 S -*O. Therefore, on substituting (5.12) and (5.13) into 
(5.14) and (5.5), we conclude that a = § , /? = | and that, within the inner region 0 ^ 
Y & 0(1), the particle concentration profile satisfies, as 0, 


2/?(<Z> 0 ) 8tf> 


$yYx 1I2 \{<p 0 ) ef 






which becomes, on integration, 


(5.15) 


e<E,_ fl/> 3(1-0.) 
dY tf> 0 (2+e 8 - 8<p «)’ 


(5.16) 


where use has been made of (5.6), (2.7) and (2.12). The solution of (5.16) is subject to 
the additional constraints that 


<Z> 0 = 0 for Y^S(X), (5.17) 

with S given by 

which follows from (5.14) plus (5.12) and (5.13). Note that since S is finite, this inner 
solution leads to a discontinuity in the particle concentration gradient at Y = S which 
can be rendered continuous by constructing an appropriate outer solution, centred 
around Y = S, to be followed by the usual matching procedure. To a first 
approximation, however, the existence of this outer solution does not affect the particle 
concentration distribution within the inner region except in the immediate neigh- 
bourhood of Y = S. 

Figure 5 shows the evolution of the particle concentration profile as given by the 
solution of (5.16) subject to (5.17) and (5.18). The similarity between this graph and 
figure 4(b) is obvious. The close agreement between the asymptotic solution as <f> s ->0 
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Figure 5. Evolution of the particle concentration profile in the region above the plate as calculated 
via the zeroth-order asymptotic solution as 0 S ->O, equations (5.17)-(5.19), at different X = X(f>* /b 
locations downstream: X = 0.0025, 0.007, 0.016, 0.03, 0.051, 0.084, 0.131, 0.195, 0.271, 0.343, 0.369. 
X cr = 0.369. 



Figure 6. Comparison of the evolution of the particle concentration along the top face of the plate 

as obtained via the zeroth-order asymptotic solution as (j> H >0 ( ), and as calculated via the exact 

solution for different values of <j> K \ , 0.2; • ■, 0.1; 0.05; , 0.02. 

and the numerical solution is further illustrated in figure 6 where (j) w , the particle 
volume fraction at the wall, is seen plotted vs. X for different values of cf) s . 

At the opposite extreme, the structure of the solution is easier to capture. 

To begin within, as seen from the evolution of the particle concentration profile with 
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Figure 7. Evolution of the particle concentration profile in the region above the plate as calculated 
via the exact solution at different locations downstream: X = 0.004, 0.008, 0.012, 0.016, 0.020, 0.024 
for <j> 8 = 0.5. X cr = 0.025 for this value of <f> s . 


X when <j) s = 0.5, shown in figure 7, the inflection point remains close to the plate even 
with increasing X. Besides, since X cr is now very small, the similarity solution presented 
in §4 applies over a larger portion of the appropriate range of X than when was 
smaller. This is seen in figure 8 which depicts the variation of <j) w with X according to 
the two solutions for the same value of <j> s = 0.5. Of course, d$ w /dX->ao as X-> X cr , 
and hence the deviations from the similarity solution become significant in this region, 
figure 8. The asymptotic behaviour in this limit can be recovered by noting that 

= + as (5.19) 


where 


X = 


X 


Y = 


Y 

(K-tsf 


(5.20) 


are the relevant stretched coordinates in this region. Further, by requiring that all three 
terms in the particle balance equation (5.5), namely convection, diffusion and 
sedimentation must be retained as <j> s -* (j ) m , we obtain that a = §,/? = §, as well as the 
governing equation for *2: 


1 0*2 


X 1/2 dX 


2/p 

HJo 


*2 2 dn 


H — 


1 3*2 


) 2 X 3l *dY 


l&M - >¥§H 


*2 2 0 2 *2 2*2 /0*2\ 2 

X m dY 2 X \0 J7 


- 2 * 2 ^?, 

0F 


(5.21) 


with 


°YL = -X^ at F = 0, 

0F 

*2 > 1 as 7-> + oo, *2->l as X^-0. 


(5.22) 


and 
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Figure 8. Comparison of the evolution of the particle concentration along the top face of the 
plate as obtained via similarity and as calculated via the exact equations for = 0.5. 


t. 

T rr 

& 




0.01 

15.000 

0.15 

1.485 

0.50 

0.0248 

0.02 

9.120 

0.20 

0.935 

0.53 

0.0125 

0.05 

4.497 

0.25 

0.580 

0.55 

0.0060 

0.08 

3.000 

0.30 

0.350 

0.57 

0.0014 

0.10 

2.415 

0.40 

0.114 

0.575 

0.0006 


Table 1. Length along the top face of the plate needed for <j> w to reach its maximum possible 

value, <j> m = 0.58, as a function of (j ) s . 


The variation of X cr is depicted in a log-log scale in figure 9 while the numerically 
computed values of X cr are given in table 1. <fi s /((f> m — <p s ) is used as the independent 
variable in figure 9 in order to facilitate comparison with the asymptotic results, 
X cr — O.36907 4/5 as ^ s ->0 and X cr — 0.352 (0 m — as (j> s ^<l> m > The factor 0.369 
was obtained by solving (5.1 6)— (5. 1 8) and corresponds to the value of X for which 
0 O (T = 0) becomes whereas the factor 0.352 corresponds to the value of X for 
which Q(Y = 0) vanishes and is obtained by fitting the asymptotic result (5.20) in 
the numerically computed values for X cr , given in table 1, in the limit A 

parameter of some practical interest is the dimensional critical distance from the leading 
edge, * cr , at which (f> w = <f> m . In view of (4.13), we have that 



where 


5 = 




Y 11 

1 




(5.24) 


which is a function of only and of the relative density ratio A p/p f . 
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Figure 9. Variation, in a log-log scale, of the distance from the leading edge along the top face of 
the plate at which <j) w = <j> m as a function of The continuous curve is produced by 

connecting the numerically obtained values whereas the two broken curves are produced by the 
asymptotic results as <f> 9 ^0 ( ), and as <f> K ^<f> m (---). 



Figure 10. Variation of 5, as given by (5.24), with <j> g for A p/p f = 0.2 ( ), and 1.0 ( ). 

The product of the last two terms in (5.24) is a montonically decreasing function of 
<fi s and approaches a constant value when as expected on account of the 

asymptotic behaviour of the exact equations in this limit derived above. On the other 
hand, since the first term in (5.24) is a rapidly increasing function of <j > s , the complete 
expression for S is a monotonically increasing function of beyond <fi s ~ 0.1, even 
up to cj) s = <f> m9 figure 10. This implies that when the suspension is densely packed, 
shear-induced diffusion is sufficiently dominant to prevent the particles from 
accumulating along the plate even for large distances from the leading edge. It should 
also be noted however that, owing to the difficulty of constructing reliable numerical 
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Figure 1 1. Evolution of the particle concentration profile in the region below the plate, as obtained 
via similarity, ( b ), and as calculated via the exact equations, ( a ), at different locations downstream: 
X = 0.008, 0.014, 0.02, 0.026, 0.032 for = 0.2. T cr = 0.032 for this value of 0,. 


solutions for (j) s > 0.575, the curves in figure 10 were generated in this interval by 
multiplying the first two terms of (5.24) with the asymptotic form of the product 
X cr ^(0s) 6/5 ~ O.352(0 m — ^) 6/5 A(^ g ) 6/5 . Similarly, since as was explained earlier in this 
section reliable numerical solutions could not be generated for <f) s below 0.01, the 
asymptotic expression X cr ~ 0.369^“ 4/5 was used in the interval 0 < <fi s < 0.01. The 
logarithmic scale used in figure 10 provides a better resolution of the curves near the 
endpoints of the interval 0 ^ <j> 8 ^ $ m . Needless to add, as the variation of S 

with <j> s is very sensitive to the choice of the parameter functions /?,/ and A used in our 
model. In fact, since the asymptotic forms of these parameters are at present not known 
with any degree of confidence, the asymptotic results given above are, of course, subject 
to revision as more reliable expressions for /?,/ and A become available for high particle 
concentrations. 

5.2. Results and discussion for the region below the plate 

Below the plate and attached to it there is a region, termed the depletion layer, where 
the particle concentration continuously increases from its value at the plate, (j ) w , to that 
in the free stream, <j > s , owing to the continuous outflux of particles sedimenting towards 
the bulk of the suspension. As a result, the particle volume fraction along the underside 
of the plate, decreases from <j> s at the leading edge until total depletion is achieved. 
Not surprisingly, near the leading edge, the concentration profiles as predicted by the 
similarity solution and those calculated via the numerical solution of equations 
(5.5)— (5.6) are in very close agreement. However, as X increases, the particle 
concentration, as obtained from the exact numerical solution, decreases rapidly with 
increasing distance from the wall and quickly reaches the constant value, <j> s9 at the edge 
of the layer, i.e. the thickness of the depletion layer is smaller than predicted by 
similarity; this behaviour is shown in figure 11 when <j> s = 0.2. In addition, as shown 
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Figure 12 . Comparison of the evolution of the particle concentration along the bottom face of 
the plate as obtained via similarity and as calculated via the exact equations for <f> s — 0.2. 



Figure 1 3. Variation, as a function of <p g , of the distance from the leading edge along the bottom face 
of the plate at which (f> w first vanishes. The continuous curve is produced by connecting the 
numerically calculated values whereas the broken curve corresponds to the asymptotic expression as 
0 .^ 0 . 

in figure 12, the particle volume fraction at the wall, <f> w9 decreases with X much more 
rapidly than predicted by the similarity solution. It should also be noted that the 
concentration profiles no longer have an inflection point and that, as <f> w ^ 0, very large 
normal gradients are observed near the wall accompanied by an unbounded growth of 
the longitudinal gradients; see figures 11 and 12. 

This behaviour can again be explained by an argument similar to that used earlier 
to account for the flow situation above the plate. Specifically, as in that case, the 
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K 

& 

X" 

0.0100 

0.0013 

0.1000 

0.0180 

0.0125 

0.0017 

0.1500 

0.0270 

0.0250 

0.0038 

0.2000 

0.0340 

0.0500 

0.0090 

0.3000 

0.0450 


Table 2. Length along the bottom face of the plate needed for <j) w to vanish, as a function of <f) g . 


variation of the diffusion coefficient which is now lower than its constant value, /?($5 S ), 
used in the similarity formulation, is the primary factor which determines the evolution 
of the particle volume fraction along the underside of the plate, hence the larger 
particle depletion rate observed in figure 12. Finally, according to (5.10), the second 
derivative of the particle concentration is negative along the plate, and since it remains 
negative throughout the layer this excludes the possibility of having an inflection point. 
Thus the situation is opposite to that in the region above the plate, where the inflection 
point moves away from the plate, and as a result the thickness of the depletion layer 
is smaller than predicted by similarity. 

The distance from the leading edge, X cn at which the concentration at the wall first 
vanishes is plotted in figure 13 as a function of <f> s with the actual numerical values given 
in table 2. A logarithmic scale is chosen to represent the data so that the comparison 
with the asymptotic result as ^ s ->0 becomes clearer. Again an asymptotic analysis that 
is uniformly valid within the depletion layer can be constructed as and by 

requiring that convection should balance shear-induced diffusion and sedimentation 
it can easily be shown that, in this limit, an 0(<^ /15 ) boundary layer forms in the 
F-direction below the plate, that extends an 0(<$j /5 ) distance from the leading edge. 
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